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ABSTRACT 



ON 
ON 
0> I Already slightly eccentric orbits, such as those occupied by many old stars in the Galactic 

disk, are not well approximated by Lindblad's epicycle theory. Here, alternative approxi- 
mations for fiat orbits in axisymmetric stellar systems are derived and compared to results 
from numeric integrations. All of these approximations are more accurate than Lindblad's 
classical theory. I also present approximate, but canonical, maps from ordinary phase- 
space coordinates to a set of action-angle variables. 

Unfortunately, the most accurate orbit approximation leads to non-analytical R{t). 
However, from this approximation simple and yet very accurate estimates can be derived 
for the peri- and apo-centers, frequencies, and actions integrals of galactic orbits, even 



VO ' for high eccentricities. Moreover, further approximating this approximation allows for 



an analytical R{t) and still an accurate approximation to galactic orbits, even with high 



Q>^ . eccentricities. 
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1. Introduction 

The knowledge of the stellar orbits is a condition for de- 
tailed understanding and interpretation of the dynamics 
of stellar systems such as globular clusters and galaxies. 
Unfortunately, even for fairly simple models for the over- 
all mass distribution in such systems, the stellar orbits 
are not accessible in a closed functional form, and numer- 
ical integration of the equations of motion is necessary. 
With the progress of computer technology this ceases to 
be a big problem. Nonetheless, it is desirable to be able 
to approximate stellar orbits analytically, in particular it 
is often useful to have simple and yet accurate estimates 
for the peri- and apo-centers of the stellar orbits and for 
their frequencies and action integrals. 

This is often done by first considering a symmetrized 
version (spherical or axisymmetric) of the stellar system 
in question, and then approximating the orbits by per- 
turbation theory. Orbits in axisymmetric galaxies, con- 
serve a component L of the angular momentum as well as 
the energy E. Orbits which are restricted to the equato- 
rial plane are always regular. Their equation of motion 
can be integrated by first obtaining the radial motion 
R{t) from 



t = 



PR 



dR' 



2[E-^{R')] -L^/R'^ 



and then the azimuthal motion by 

'* 'dH\ .., /■* Ldt' 

'dL 



m 



dt' 



(1) 



(2) 



R{t') J ^ (*') 

where $(i?) denotes the gravitational potential, while 
H{R, L,pji) is the Hamiltonian. Unfortunately, the inte- 
gral (yj) can be solved analytically only for very a few po- 
tentials, notably that of a point mass, and the harmonic 
and the isochrone potentials, neither of which gives a 
good description of the potential of real galaxies. There- 
fore, one often employs an approximate solution to (El), 
namely Bertil Lindblad's (1926) epicycle theory, well- 
known from the textbooks (cf . [Binney fc Trcmaine 1987| ) . 
This approximation may be derived from equation (111) by 
expanding the argument of the square root in the denom- 
inator into a quadratic in R, yielding an explicitly soluble 
integral. However, as we shall see below, the argument of 
the square root is not well approximated by a quadratic 
in R, and neglect of the higher-order terms causes the ap- 
proximation to be useful only for orbits with quite small 
eccentricities. 

There are two different routes to improve on the clas- 
sical epicycle motion. The first one is an extension to 
higher-order perturbation theory, which can be done in 



various ways (cf. Lindblad 1958; Kalnajs 1979). The al- 
ternative, which I will pursue in this paper, is to stay at 
first-order perturbation theory but use a better approx- 



imation for the integrand in (||) . A recipe for doing so is 
presented in Section 0, while some applications, i.e. con- 
crete orbit approximations, are given in Section S. Lind- 
blad's and Kalnajs' (1979) epicycle theories are recovered 
as the two simplest cases. There are three other approx- 
imations of interest, all of which are more accurate than 
these epicycle approximations for the important case of 
a flat rotation curve. 

I also consider the approximate integration of the az- 
imuthal motion (0) with particular emphasis on obtain- 
ing a canonical map to the azimuthal angle variable 0^. 
This allows for a approximate but canonical map to a 
set of action-angle variables (J^, L, 6*^, 6*0) and guaran- 
tees that the resulting approximated phase-space flow if 
incompressible. Such an approximate mapping is most 
simple for the classical epicycle theory (though, to my 
knowledge, not given in the literature). 

In Section y, a new explicit orbit approximation is 
given, which has been derived from the most accurate, 
but unfortunately implicit, approximation of Section ^. 
As comparisons with numerically integrated orbits in 
Section ra show, this approximation is indeed very ac- 
curate, even for highly eccentric orbits. The paper is 
summarized and concluded in Section ^. Appendices A 
and B give useful relations for circular orbits and power- 
law models, respectively. 

2. A Formalism for Approximating 
Stellar Orbits 

In this section, a general recipe is presented for approx- 
imately solving equation IM) for R{t) and equation (0) 
for (t>{t). I also derive a way for obtaining an approxi- 
mate map between ordinary phase-space coordinates and 
action-angle variables. Worked examples, i.e. concrete 
orbit approximations are presented in Section 0. 

2.1. The Radial Motion 

2.1.1. Modifying the Time Integral 

The key idea in obtaining an approximation superior to 
classical epicycle theory is to transform the integral (|^) 
into a form which is better adapted for replacing the 
argument of the square root by a quadratic. Consider 
re-writing equation (n^) as follows 

dR _ r ^{R) dR f a;" dx 

PR 



t^ 



(3) 



aR)PR J aR)PR' 

where £,{R) is a positive definite function of R. The re- 
lation for x{R) follows upon integrating £,{R) dR — a;" dx 
with integer n. We now introduce the auxiliary variable 
77 by 

^=^\ (4) 

d?7 
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which essentially amounts to a new time variable. For 
n > 0, equation (0) means that, for fixed drj, the time- 
step dt is smaller near peri-center than near apo-center. 
With (^), equation (||) becomes 

dx 

(5a) 



V = 



Vn^ 



with 



E-'^{R) 



2i?2 



(5b) 



The new integral ( pa| ) is formally identical to the orig- 
inal (El), and we might derive an approximate solution 
for xlrj) in the same way Lindblad's epicycle theory is 
derived from (m. Approximating 

Y^X^-a^{x-xof, (6) 

obtained, for example, by Taylor expansion, we find 

X ^ xq{1 — e cos arj), e — X/{axQ). (7) 

With this form of x(ry), equation (^ is easily solved for 
t{ri). However, ri{t) cannot be obtained analytically un- 
less n — 0. 

2.1.2. What Functions ^{R) are Useful? 

Two choices, ^=1 and S,= R, are of particular interest. 
For the first, 



Y ^2E- [2^{R) + L^/R^]. 



(8) 



For orbits with the same angular momentum L, the Tay- 
lor expansions of Y are identical apart from the additive 
constant 2E. The maximum of (pi) occurs at the radius 
i?i = Rc{L) of the circular orbit with angular momentum 
L. At this radius Y = 2J\E with 



AE = E-E,{Rl), 



(9) 



where Ec{R) is the energy of the circular orbit at radius 
R (see equation A3 for a definition). 
For the second choice, £, = R, 



Y = 2i?2 [E - $(i?)] - L^ 



(10) 



This time the Taylor expansions of Y are identical, apart 
from the additive constant — L^, for all orbits with the 
same energy E. The maximum of (M) occurs at the 
radius Re = Rc^E) of the circular orbit with energy E. 
At this radius 



Y ^ AL^ = LKRe) - L^ 



(11) 



where Lc{R) is the angular momentum of the circular 
orbit at radius R (see equation A2 for a definition). 

For all other choices of £,{R), the radius at which Y 
becomes maximal depends on both E and L in an essen- 
tially non-trivial way. 



2.1.3. The Modified Equation of Motion 

We might also give the equation of motion that corre- 
sponds to the approximation introduced in Section |2.1.l| . 
The integral ( |5a| ) solves 



1 dY 



(12) 



-— = -a [x ~ xo) 



drj^ 2 dx ' 

while the approximation (|7|) actually solves 

d^x 
d?7 

with X being a constant of integration. Note that equa- 
tion (n2h in conjunction with equation (0) yields 
.. _ 1 9(ra;-2« 
^~ 2 dx 

which is equivalent to the familiar R = —OH/dR. The 
approximated equation of motion dlS) together with 
equation (Q) is equivalent to i? == — (dH/dR) with the 
approximate Hamiltonian 

1 



(13) 



(14) 



H = E 



Pi 



2^2 



[X^-a'{x-x^f]+g{E,L). (15) 



Note that in general H depends on E and L as param- 
eters of X, xo, a, and, possibly, x{R) and £,{R). I have 
added a term g(E,L), where g is an arbitrary function, 
because this addition does not change the dynamics for 
Rit). 

2.2. The Azimuthal Motion and 
Action-Angle Variables 

After obtaining an expression for R{t), the radial an- 
gle, 9]i, is easily identified as the term, linear in time, 
of which R{t) is a 27r-periodic function. With the above 
formalism, the radial action is 



^2 /•^'^ sin^ary x" dry 



^'^ ^dRV fdxV dr] 



dx 



d-qj dt 



dr] 



27r 



/o e{R{^)) ' ^^^^ 

where x = x{rj) is given in equation (R). For this integral 
to be soluble, £,{R) must not be too complicated when 
expressed as function of x. 

The azimuthal motion may be obtained by inserting 
the approximation for R{t) into equation (g) and inte- 
grate, if necessary, after further approximating the in- 
tegrand. However, the approximate orbits generated in 
this way do not necessarily produce an incompressible 
flow in phase-space, as Liouville's theorem demands. In 
other words, the part of (j){t) which is linear in time is 
not, in general, the angle 0^ canonical conjugate to the 
angular momentum L at fixedn {Jb.,0b.). When using 



Note that at fixed (i?, Pij), the azimuth itself is conjugate to L. 
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such an approximation, spurious effects may result from 
this lack. Note that in particular the textbook variant of 



Lindblad's epicycle theory (cf. Binncy & Tremaine 1987 
p. 124) suffers from this problem; see §3.1.1 below for a 



variant of Lindblad's theory that generates a Hamilto- 
nian and hence incompressible phase-space flow. 

An incompressible phase-space flow or, equivalently, a 
canonical map to action-angle variables {Jr, L, 9r, 9^), is 
best generated from the approximate Haniiltonian (Eq). 
That is, the azimuthal motion is approximated by inte- 
grating (p = dH/dL, and the angle 9^ is identified as the 
contribution to (j}{t) which is linear in time. 

A complication arises, because the function g{E,L), 
which appears as additive term in H, does affect (/> by 
adding the constant dg/dL to the azimuthal frequency, 
even though it makes no difference for R{t). Thus, for 
the sake of obtaining the relation 6*0 (t) (but note the 
canonical map itself), we must specify the function g. 
This can be done by requiring that 



dH 



= ^R, 



(17) 



where lor is the radial frequency emerging from the ap- 
proximation for R{t). Inserting the approximation (|^) 
into ([l5|), one finds that H — E + g. Together with Jr, 
obtained from equation (iq), the requirement (nW gives 

fdJR^ 






i^R 



\ dE 



1. 



(18) 



With this choice of g{E,L), the approximated orbital 
frequencies satisfy the identity 
du}R\ d'^H f duj^\ 



dL 



Jr 



dJRdL ydJRj^ 



3. Worked Examples 

3.1. Recovering Known Approximations 

3.1.1. Classical Epicycle Theory 

Lindblad's (1926) classical epicycle theory is recovered 
from the above recipe for ^ = 1 and n = {]. The Taylor- 
expansion of Y around R^ reads, with 5r = R — Rl, 



Y = 2A^ -K^^R- 






3k^ 

~r" 



Rl 



^+0(4), (20a) 



where the epicycle frequency k, defined by equation (A5), 
is evaluated at R^, and hence is a function of angular 
momentum. For the important case of a fiat rotation 
curve. Figure 1 shows Y__2E and its approximation by 
the quadratic part of (20a). The large error of this ap- 
proximation is related to the neglected third-order term 
in ( ^Oa| ), which for stellar systems is never small, since 
always d/t^/di? < 0. The radial motion is approximated 

by 

R{t) = Rl {I -e cos 0r), (20b) 



(~1, n —0: classical epicycle theory 



-SE^{L) 



Y-2E=-2'I>(R)-1A 
quadr. approx. 




R,,(;-) x=R 

Fig. 1. — Lindblad's classical epicycle theory for the logarithmic 
potential (flat rotation curve). Th e solid line gives Y{R) (pi) while 
its quadratic approximation (20a) is shown as broken line. The 
thin horizontal line corresponds to the energy of an orbit with 
L = 0.9Lc{E). Note the strong asymmetry of Y{R). 



e^V2AE/{RLn), 
Jji^^KRle^ ^ AE/k, 



(20c) 
(20d) 



where 9r = nt. The approximate Hamiltonian is, with 
Ec^Ec{Rl), 



H = Er 



[PR 



,2r2- 



Ec + kJ^ 



R, 



(20e) 



specifically, g = 0. The azimuthal motion follows from 
integrating — dH{R,pR, L)/dL: 



e^-^Rlsm29R 
dL 



(20f) 
(20g) 



(p = 9^ + e-f sm 9r -i 

UJ4, = n{L) + -tj-Jr, 

where 7 = 2J7/k is evaluated at Rl , while 9^ = uj^t. It ap- 
pears that these formulae (pOl) and (20g) for the canon- 
ical map (f>{9^) and the corresponding azimuthal fre- 
quency io^ that are consistent with the radial motion 
(20b) are not in the literature. Note that the error of 
equations (20i) and (|20g|) is still 0{e'^), even though they 



contain terms of order e (x Jr. Omitting those terms 
gives the textbook result for (/)(i), but does not yield a 
canonical map to the angle 9^ and results in a compress- 
ible phase-space flow. The azimuthal frequency (20g) 
and the radial frequency ujr = k{Rl), satisfy the idcn- 
tity (|l|). 

3.1.2. Kalnajs' Epicycle Theory 

In a little-known four-page paper, Kalnajs (1979) intro- 
duced an improved epicycle theory, which can be derived 
from our recipe with ^=ii? and n = 0, giving x — R^. 
The Taylor expansion of Y (equation |lO|) around its max- 
imum is, with Sx = R^ ■ 



-R|> 



Y = AL' 



1 



dK 



2 1 



di?2 



6^ 
12 



0{5t). 



(21a) 



Here, k and its derivative are evaluated at Re, i.e. they 
are functions of energy. Y {x) 4 - L^ and its approximation 
by the quadratic part of (21a) are shown in Figure g for 
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Lim 



i^R/2. n=0: 



Kalnajs' epicycle theory 

Y+L^ = 2R^[E--t>(R)] 

quadr. approx. 




RfCB) 



R^...(P) 



Fig. 2. — Kalnajs' improved ppinynle theory for the logarithmic 
potential. T he s olid line gives Y{x) (110) while its quadratic ap- 
proximation ( piaj ) is shown as broken line. The thin horizontal line 
corresponds to an orbit with L = 0.9Lc{E). 

the logarithmic potential. The resulting approximation 
for R{t) is 

R{t) = Re {I -ecoseny/'^, (21b) 

e^jy/l-L^/Ll, (21c) 



Jr^ 



27 



1- VT^ 



(21d) 



where 9f{ = Kt, j — "/{Re), and Lc = Lc{Re)- As pointed 
out by Kalnajs, the neglected third-order coefficient in 
(21a) is always smaller than that in (EOq), implying that 
his theory is generally more accurate than Lindblad's. 
This is also evident when comparing Figures 1 and 0, 
where the horizontal lines correspond to the same orbit. 
Kalnajs also remarked that his theory becomes exact for 
a harmonic potential, which arises from a uniform mass 
distribution and, at the time he wrote the paper, was 
thought to be a good description for the inner parts of 
galaxies. 

Integrating 4> — dH/dL exactly gives 



(21e) 




(21f) 
(21g) 



with 9^ —Lo^t and fl — fl(RE)- For the harmonic poten- 
tial, ^ — ^oj^R^, K — 2uj and g — 0, while for a flat rota- 
tion curve, Vc = vo, 



dg_ 
dL 



1 



1 



x/T— 



(21h) 



2i 
In equation fple 
A{x,y) = arctan (a; tany) — arctan(tanj/) (22) 

is a continuous, periodic function of y with period tt. 



ms) 



^=R, n — \: Isoclrirone approxirnalion 

y+L2 = 2R2[£'-*(fi)] 

quadr. approx. 




.„.„(E) 



Fig. 3. — The isochrone approx ima tion for the logarithmic po- 
tential. The solid line gives Y{x) ( |lO[ ) while its quadratic approx- 
imation (E3) is shown as broken line. The thin horizontal line 
corresponds to an orbit with L = 0.9Lc(E). 

3.2. Approximations with ^ = i?, n = 1 

3.2.1. The Isochrone Approximation 

For this choice of ^ and n, one finds x = \/R?+l? where 

The Taylor expansion of 
with 



h is an arbitrary constant. 

Y{x) around its maximum at xo = \/Re + ^^ is 

Sx =x ~ xq and k — k{Re), 



Y^AL' 



.2 ,.2 c2 



2x1 dK^ 



xo^n' + ^^Sl+OiSi),{23) 

while the resulting expressions for R{t) and (p{t) can be 
found in Appendix O The neglected third-order term in 
(p3) can be made to vanish identically by the choice 



-Ri 



din/ 



21-1 



dlni?2 



1 



(24) 



For all realistic stellar systems, except the harmonic po- 
tential, for which Kalnajs' theory is exact, equation (p3) 
results in well-defined b^ > 0. Figure actually plots 
Y{x) + L^ for a flat rotation curve and b chosen according 
to (|^ . A comparison with Figure g shows the qualita- 
tive improvement over Kalnajs' epicycle theory. This ap- 
proximation becomes exact for the potential for which b 
is constant, resulting in Hcnon's (1959) isochrone sphere 

""'' (25) 



$ 



b + VWT¥' 



3.2.2. 



The Keplerian Approximation 

I will also consider the simple case 6 = 0. The resulting 
approximation is exact for orbits in the potential of a 
central point mass, corresponding to (|2^) for b — Q. For 
a flat rotation curve. Figure ^ plots Y{x) + L"^ and its 
quadratic approximation. A comparison with Figure |^ 
shows that for this case the Keplerian approximation is 
better than Kalnajs' theory. The relations for the radial 
motion are familiar from celestial mechanics (for which 
7 = 2): 

R{t) ^ Re {I - e COB rj), (26a) 

6]l — Kt — rj — e sin rj, (26b) 



Walter Dehnen 



Lm 



(=R, n = l, b-0: 



Kepleriari appi'oxirnation 

Y+L^ = 2R'^[E--t>(R)] 

quadr. approx. 




njE) 



n^.J'^) 



Fig. 4. — The Keplerian approximation for the logarithmic po- 
tential. The solid line gives Y{x) (|lO|) while its quadratic approx- 
imation (bsl) is shown as broken line. The thin horizontal line 
corresponos to an orbit with L = 0.9Lc{E). 



2Lr 



Jr 



1 



1-v^ 



(26c) 
(26d) 



The relation between azimuth and time that results from 
integrating ((> = dH /dL is, with 6*^ — oj^t, 



72 - 4e2 



esiuT; 




uj^ — Vl 



1 — 46^7 2 Qg 



l-e2 



AE 



EJRl) 



dL' 
dluK" 



dlni?2 



1 



VT^^ 



,(26e) 

(26f) 
(26g) 



For the potential generated by a central point mass, 
g — 0, while for a flat rotation curve. 



dg 
dL 



L 



1 - 



1 



VT 



(26h) 



3.3. ^ oc R and n — j(RE) -1 

The approximation that results from our recipe for this 
choice of ^ and n, giving 

X = {R/REf'\ (27a) 

is of particular interest. For the harmonic and Kepler 
potential, this yields, respectively, Kalnajs' epicycle the- 
ory and the Keplerian approximation, i.e. the exact so- 
lutions. For other potentials, n is non-integer, and con- 
sequently tirf) is not analytical. However, one may still 
give R and pn as functions of -q and evaluate Jr. Here 
are the basic results (with 5^=x—l and Lc — Lc{Re)) 



Y ^ 


-Ll 


^2 x2 ( „, 1 1 7 dliiK^ 
[^ ^x [1 J^+3dlnfl2 


e = 


= Vi 


-L2/L2, 


R = 


= i?B(l-ecosr;)''/', 


Kt - 


=/;'( 


1 — ecosTy') dr;'. 



'0{5t 



%{Lc 



(27b) 

(27c) 

(27d) 

(27e) 
(27f) 



(<^-R , ri=y- 1 



Lm 



Y+L^ = 2R^[E-<S>{R)] 
quadr. approx. 




x^{R/Rj'/y 



...m 



Fig. 5. — The approximation arising from ^oc R and w = 7 — 1 for 
the logarithmic potential. The solid line gives Y(x) (|lO[), while its 
quadratic approximation {27b) is shown as broken line. The thin 
horizontal line corresponds to an orbit with L = 0.9Lc{E). 



Note that the neglected third-order term in (27b) is usu- 
ally smallR: in the case of power-law rotation curves, 
VcOcR^, its coefficient becomes maximal in magnitude 
(at —0.0572 . . .) for f3~0. For this worst case of a flat 
rotation curve. Figure plots Y{x) and its quadratic ap- 
proximation. Evidently, this approximation is superb, 
even for small angular momenta. In particular, it always 
gives a peri-center radius of zero for L — 0, which all the 
other approximations fail to do (apart from the cases 
where they are exact). Note that the radial frequency 



Wfi 



,7-1 



(|L|/L,)^-^P,_i(Le/|L|) 



(27g) 



where P,y(x) denotes the associated Legendre function, 
is larger than k, unless L — L,., 7=1, or 7 = 2. 
Integrating (j> = dH /dL gives 



= 6*0 + 7 


2 


""+.4 


( /1 + e v\ 
VVl-e'2| 


; 


(27h) 


..^o^ + ||, ,2n, 


JE^iRr.) ^ « 


^ ' dhlfl.2 l^^ La) 


- 1) •(27j) 


with 6*0 = iv^t and Or — urI. For a flat rotation curve, 


dg n 


'^R / 
. K 


(27k) 


dL Vl- 


-e2 



4. A New Orbit Approximation 

The remarkably accurate orbit approximation presented 



in 53.3 above has the unfortunate drawback, that the 



time t cannot explicitly solved for, because the integral 



in equation (27e) is generally not analytically soluble. 



However, one may approximately solve this integral, by 
replacing its integrand by a finite powers series in e. If 



Replacing (27a) with x"' = (iJ^-f 6'^)/i?|, and choosing the constant 



b such as to set the third-order term to zero, as in the isochrone 
approximation, does not work, since it requires b^ < 0. 
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i;{E) 



a new approxirnaiion 



y+L2-2/?^[f;-$(R)] 
Isl order 
2nd order 




x-(R/R^r/y 



,„(f^) 



Fig. 6. — The approximations described in §W for the logarithmic 
potentiah The solid line gives Y{x) (|lO[), wh ile th e dashed and dot- 
ted lines correspond to the approximation ( 28a ) for fcmax = 1 and 
^max = 2, respectively, for the two orbits with L/Lc{E) = OA and 
0.9 (for the latter the approximations can hardly be distinguished 
from the exact relation). 

the highest power retained in this procedure is e'^™'"' , this 
technique amounts to approximating, in the forniahsni of 



Y = R'pj, « Li (e^ 



-\ 2 



V k J 



bl 



(28a) 



E^max 
fc = 

For the case of a flat rotation curve. Figure H plots Y{x) 
and the approximation ( |28a| ) with fcmax = 1 (dashed) and 
fcmax = 2 (dotted) for the two orbits with LiLc = 0.4 and 
0.9. Evidently, these approximati ons a re not much worse 
than their parent, i.e. equation ( 27b ), and at least as 
good as the isochrone approximation above. The radial 
motion, R{t), is given by 

R = RE{l-ecosTjy^'^, (28b) 

e^^yi-L^L-^, (28c) 

Kt ■ 



rj — e (7— 1) sin 77 



+ ie2(7-l)(7-2)(77 + isin277) + 
LUR^K [l + ie2(7-l)(7- 2) + ...]"'. 



(28d) 
(28e) 

Equation (28e) implies ujr>k since always 1<7<2, 
and is, for power-law rotation curves, even more accurate 
than its parent ( |27^ , i.e. including terms 0{e^) would 
degrade the accuracy of the new approximation. Note 
that the corresponding radial action cannot be given in 
closed form. The azi muth al motion may be approxi- 
mated using equation (27h) and replacing dg/dL by its 
result for a flat rotation curve (which for power-law mod- 
els makes an error of the order e*^): 



(f) = w^i -I- 7 



W0 




(28f) 



(28g) 



The relation (28d) for t{r]) is generally not invertible. As 
a consequence, one must either (i) use rj rather than t 



as independent variable, (ii) solve (28d) numerically for 
77(f), or (iii) use the approximate inversion 

T] w LURt + e(7 - 1) sinujRt 

+ le\j-l){3j-2)sin2ujRt. (29) 

Note, however, that using this relation is not equivalent 
to using (28d), rather it generates a different orbit ap- 
proximation. 

As we will see below, the orbit approximation defined 
by equations (Eq) is remarkably accurate down to high 
eccentricities. However, it is not Hamiltonian, i.e. the 
phase-space flow generated is in general not incompress- 
ible. It deviates from this ideal by an amount of the 
order 6*^™^'=+-'^. 

5. Comparison with Exact Orbits 

We now compare the various orbit approximations with 
numerically integrated orbits in the logarithmic poten- 
tial, $ = Uq ln(_R/_Ro), which supports a flat rotation 
curve. This potential is scale invariant, i.e. all orbits 
with the same ratio \L\/Lc{E) but different energies are 
identical, apart from a scaling relation or the sense of 
rotation. 

5.1. Testing the Radial Motion 

For eight orbits with \L\/Lc{E) between 0.5 and 1, Fig- 
ure compares the plots of R vs. pr (surfaces of sec- 
tion) predicted from the various orbit approximations 
with the exact relations. The agreement in this plot 
is indicative of the accuracy of the relation for R{t). 
For a quantitative comparison with the orbit of a star 
with velocity {U, V, W) in the solar neighborhood, one 
should compare pr with the radial velocity U. The az- 
imuthal velocity V w.r.t. the local standard of rest can 
be evaluated from the ratios R/Rc{E) and L/Lc{E) to 
be V = {[L/ Lc{E)]/[R/ Rc{E)] - l)wo with vq denoting 
the local circular speed. 

Evidently, Lindblad's classical epicycle theory be- 
comes significantly inaccurate for |L|^0.98Lc, corre- 
sponding in the solar neighborhood to stars with |C/| ^ 
45kms"i or \V\ > 32 kms^^ (with wq ~ 220kms"i). A 
significant fraction of old-disk stars in the solar neighbor- 
hood have velocities in excess of these values (cf. Dehnen, 
1998), i.e. the classical epicycle theory cannot be safely 
used in quantitative studies of the local Galactic disk. 

Kalnajs' epicycle theory is applicable down to \L\ « 
0.95Lc(£')- Theories with n = l are better than the 
epicycle theories (n = 0). The Keplerian approximation 
predicts too large peri-centers for |L| ^0.85Lc(-E'), and 
is thus just capable of describing R{t) for stars in the 
high-velocity tail of the local stellar disk distribution. 
The isochrone approximation is superior to these, giving 
accurate approximations down to \L\ fnO.TLdE). 
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exact orbit 

classical epicycle theory" 
Kalnajs' epicycle theory 
Keplerian approximation 
Isochrone approximation- 
i-R.n-y{E)-l 
equations (28), 1^'^ order 
equations (28), 2"'^ order 




Fig. 7. — Surfaces of the R-pR section of phase space: comparison of exact orbits in the logarithmic potential, which supports a flat 
rotation curve, with the predictions of the various orbit approximation of §N and §y. The orbital angular momenta are, from inside out, 
0.98, 0.95, 0.9, 0.85, 0.8, 0.7, 0.6, and 0.5 of Lc(E). 



0.5 



exact orbit 

classical epicycle: 1st order 
classical epicycle: AAV 
Kalnajs' theory: AAV 



Keplerian approx: AAV 
Isochrone approx: AAV 
equations (23), k^^^-] 
equations (23), ^^^^3^=2 




R / R„(S') 

Fig. 8. — Surfaces of the R-{if> — ui^t) section of phase space: comparison of exact orbits in the logarithmic potential, which supports 
a flat rotation curve, with the approximations for the azimuthal motion. The orbits have the same angular momenta as in Figure M. 
Approximations which allow for a canonical map to a set of action-angle variables are indicated as 'AAV. Note that multiplying the 
j/-axis with R, these essentially are the epicycles which the orbits perform around their guiding centers. 



Also shown are the relations arising from the approx- 
imations in §3.3 and §0. Evidently, these are very good: 



the approximation derived in §3.3 and its second-order 
expansion (equations Eq with /cmax = 2) are much better 
than any of the other approximations. Even the first- 
order expansion (fcmax = 1) is comparable in accuracy to 



the isochrone approximation, which itself is second-order 
accurate. 

5.2. Testing the Azimuthal Motion 

To test the quality of the approximations for the az- 
imuthal motion, I compare in Figure H the surfaces of 
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Fig. 9. — Radial action Jji (left column) and frequency ujn (middle column) and azimuthal frequency u)^ (right column): comparison 
with exact orbits for the potentials supporting the rotation curves VcOcRl^ with /3 = — 0.2, 0, and 0.2 (from top to bottom). Note tha ; 
'jJR = f^(E) for all these approximations except Lindblad's epicycle theory and the approximation of §W with femax = 2 (equation E8e, 
second order). 



R-{(l) — Lo^t) sections of the same orbits as in Figure [^ 
with some of the approximtions. 

The approximations for which (j}(t) has been obtained 
such as to yield a canonical map to a set of action-angle 
variables, are indicated as 'AAV in the figure index. 
Among these approximation, the order of accuracy is 
similar as for the radial motion, but they are less accu- 
rate in this section of phase space than in the R-pR sec- 
tion. This is not surprising, as the approximations were 
designed to approximate R{t), and any deviation from 
the true R{t) is likely to be amplified when integrating 
for 0(t). 

Figure @ also plots three orbit approximations which 
do not yield a canonical map to action- angle variables: 
the textbook variant of the classical epicycle theory and 
the new approximation of §H for first- and second-order 
expansions. This variant of the epicycle theory is clearly 
superior in accuracy to its canonical counterpart, i.e. en- 
forcing a canonical mapping degrades the accuracy of the 
approximation. The two approximations resulting from 
equations (Eq) are much more accurate than all the oth- 
ers, though they are slightly less accurate than in the 
R-PR surface of section. Strangly, the first-order approx- 
imation is more accurate than the second-order. 



5.3. Testing Prequencies and the Radial Action 

In Figure ^, the radial action and orbital frequencies pre- 
dicted by the orbit approximations are compared with 
the exact ones for all orbits in the potentials support- 
ing nearly flat rotation curves Vc oc R^ with /? — —0.2, 0, 
and 0.2. Even for nearly eccentric orbits, the classical 
epicycle approximation under-estimates Jr and drasti- 
cally over-estimates the frequencies. The three approx- 
imations with ^ocR and integer n (Kalnajs, Keplerian, 
isochrone) have lor = k(Re), which actually is not a bad 
description of the truth. Similarly, uj^ = Q{Re) is a rea- 
sonably good estimate. The estimates for ujr from the or- 
bit approximations with canonical maps to action-angle 
variables are only useful for large L, corresponding to 
small e. 



I also plotted the estimates (28e) and (28g) for ujr and 
Lo^ respectively. While the latter becomes inaccurate at 
L <^ 0.5Lc, the estimate for cur is very accurate for all L 
(the corresponding lines are almost entirely overlaid by 
those for the exact relations). 

At Lk Lc{E)^ all orbit approximations give Jr, cor- 
rectly, insomuch that they agree in the linear term when 
expanding Jr as power series in 1 — \L\/ L^- This linear 
relation is actually the result for Jr obtained from the 
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approximation of §3.3 



6. Conclusions 

The classical epicycle theory of Lindblad (1926) is usu- 
ally derived starting from the Hamiltonian or, equiva- 
lently, from the equations of motions for a near-circular 
orbit in a spherical or flat axisymmetric potential. In 
this paper, I instead directly considered the integral (|^) 
to be solved in order to obtain R{t). The simplest way 
to approximate this integral by some closed form leads 
straight to the classical epicycle theory. In Section ^, I 
derive a general method for approximating this integral 
in a better way by manipulating the integrand before 
approximating it. 

6.1. Relation to Other Work 

Several published studies deal with improved approxi- 
mations for the stellar orbits. Most of them start with 
Lindblad's classical epicycle theory and try to improve 
it by extension to higher orders, variations of the con- 
stants involved, or both. Apart from the quoted paper 
of Kalnajs, I know of only one more study, Shu (1969), 
in which an orbit approximation is derived that does not 
incorporate Lindblad's theory. 

Shu also considers directly the integral (|l|) and ma- 
nipulates it to obtain approximate solutions. Unlike this 
work, Shu uses a (truncated) power-series expansion in 
the eccentricity, the lowest order of which yields Lind- 
blad's theory, while the next order yields an approxima- 
tion that becomes exact in the case of Keplerian motion. 
This approximation is not quite the same, though, as 
the Keplerian approximation above, for example, the ra- 
dial frequency for Shu's theory is ljii — {1 — e^)^/^K(i?L) 
with e given in equa tion ( ^Oq ), while the approximation 
proposed in Section [3.2.2| has lor = k{Re)- 

6.2. New Orbit Approximations 

From the general recipe in Section g, endlessly many 
orbit approximations may be derived. In Section 0, 
four of them are presented in addition to the classical 
epicycle theory. One has already been introduced by 
Kalnajs (1979) and becomes exact for the harmonic po- 
tential. Two other approximation yield the exact orbits 
for the potentials of, respectively, a central point mass 
and Henon's (1959) isochrone sphere, and I call them 
accordingly Keplerian and isochrone approximations. A 
comparison with exact orbits in the logarithmic poten- 
tial, which supports a flat rotation curve, shows that 
the latter two approximations are better than Kalnajs', 
which in turn is superior to the classical epicycle the- 
ory. In particular, they may safely be applied even to 
the high- velocity tail of the stellar old-disk population. 



The fourth new orbit approximation may be consid- 
ered a hybrid between Kalnajs' theory and the Keplerian 
approximation, it is very accurate, even for purely ra- 
dial orbits, for potentials with power-law rotation curves, 
VcOcR^. However, this approximation has the severe 
drawback that the time cannot explicitly be solved for. 

6.3. Action-Angle Variables 

After obtaining an approximation for R(t), I also con- 
sider the task of integrating the azimuthal motion such 
that the resulting approximation for the phase-space flow 
is incompressible, as Liouvillc's theorem demands. This 
is equivalent to obtaining an approximate but canoni- 
cal map to a set of action-angle variables {Jj^, L, 6'fj, 6^,). 
For Lindblad's and Kalnajs' epicycle theories, as well as 
for the new Keplerian and isochrone approximation, this 
map can be given explicitly. However, the azimuthal fre- 
quency LUij, that is consistent with this map cannot be 
expressed in closed form, except for the classical epicycle 
theory or otherwise for special potentials (the logarith- 
mic potential and that for which the respective approxi- 
mation is exact). 

Using an approximation for the azimuthal motion that 
does not yield an incompressible phase-space flow may 
lead to spurious effects. An example for such an approx- 
imation is the textbood variant of Lindblad's classical 



epicycle theory. In Section 3.1.1 , a version of this theory 
is presented that does not suffer from this problem. 

6.4. Estimates of Orbital Properties 

In face of the ease by which stellar orbits can be com- 
puted numerically, orbit approximations may be most 
important as a tool for quantitative estimates of stellar 
dynamical processes, in particular in dynamically cool 
stellar disks. However, a comparison of the approxima- 
tions with integrated orbits in Section showed that 
the classical epicycle theory greatly over-estimates the 
orbital frequencies and under-estimates the peri- and 
apo-center and radial action. Consequently, quantitative 
analyses based on this approximation will inevitably in- 
volve systematic errors. The size of these error may reach 
10% of more already for slightly eccentric orbits, which 
in the solar neighborhood correspond to velocities w.r.t. 
LSR in excess of about 40kms~^. 

Highly accurate estimates for the orbital quantities 
may be de rive d from the orbit approximation presented 
in Section 3.3. This approximation, even though it does 
not yield R(t) in closed form, is very accurate for orbits 
in potentials which support nearly flat rotation curves, 
like that of the Galaxy. The resulting estimates for the 
peri- and apo-centers are 



R±iE,L)^R,{E){l±e) 



7(£)/2 



(30a) 
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with eccentricity 

e^ = l-L^/Ll{E) (30b) 

and 7 = 2f}/K. Here, Rc{E), Lc{E), n{E), and k{E) 
denote the radius, angular momentum, azithmal fre- 
quency, and epicychc frequency, respectively, of the cir- 
cular orbit with the same energy E as the orbit in ques- 
tion. For orbits in potentials with power-law rotation 
curves VcO(:R^, equation ( |30a| ) under-estimates the apo- 
center by at most 1.7% (for L — 0, /3«— 0.16), while 
the peri-center is under-estimated by at most ~ 2% (for 
i«0.3Lc, /3~0.1). The corresponding estimate for the 
radial action is 



jR{E,L)^^j{E){L,iE)^\L\). 



(30c) 



For power-law rotation curves, the largest error of this 
estimate occurs on purely radial orbits {L — 0) in the log- 
arithmic potential {(3 — 0), where actually Jr — Lcy/e/2TT 
on radial orbits (equation B4), while the above approxi- 
mation predicts a value 7.5% larger. 

Another way to use this very accurate orbit approxi- 
mation is to further approximate it, in order to overcome 
its major drawback and solve for the time explicitly. This 
results in a new orbit approximation, presented in Sec- 
tion m which gives R{t) and is still very accurate to high 
eccentricities. The corresponding approximation for the 
radial frequency 



u;r{E,L) « k{E)/{1 + 4(7 - 1)(7 - 2)) 



(30d) 



is, for power-law rotation curves, accurate to better 
than 1.4% (at L = 0, /3~0.3), while the rough estimate 
ujr{E,L) « k(_E) under-estimates the radial frequency 
by at most 7.5% (at L = {] = (3). The azimuthal fre- 
quency, which may be expressed as the time average of 
LR~^, is the orbital quantity most difficult to estimate, 
because it is very sensitive to R{t) near peri-center, where 
the approximations for the radial motion are least ac- 
curate. Numerical experiments suggest that the simple 
estimate 

W0 « n{E) (30e) 

is still the most useful: it becomes inaccurate at the 
^10% level for orbits with L <^0.3Lc. 

Equations (30a) to (30c) provide accurate and yet 
simple estimates for basic orbital properties in typical 
galactic potentials. Note however, that these estimates 
are not self-consistent, i.e. they do not in general satisfy 
ujR = {dE/dJR)L or uo^ = {dE/dL)j^. 
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ful comments. This work was financially supported by 
PPARC. 



A. Circular Orbits 

This appendix summarizes the relations between the 
quantities describing circular orbits. While some of these 
relations are well-known, several others, even though 
trivial in principle, are less familiar. A circular orbit at 
radius R is defined by equilibrium between gravitation 
and centrifugal forces, resulting in the relation 

vl{R)^R-^{d^/dR) (Al) 

for its azimuthal velocity. A circular orbit is uniquely 
determined by either its radius R, angular momentum 
L, energy E, circular frequency r2, or radial frequency k, 
also known as epicycle frequency. (Radial oscillations are 
not excited for circular orbits, but their radial frequency 
is well-defined nonetheless). The definitions and some 
relations between these are 
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B. Power-Law Potentials 

Quite often one is dealing with simple power-law models 
with circular speed 



VciR) = ""0 ^ 



R 

Ro 



/3 



and gravitational potential 

2/3 



$(i?) 



2/3 \RoJ 






for /? 7^ 



for /3 = 0. 



(Bl) 



(B2) 



Here, Ro and vo are scale radius and scale velocity, and 
shall be set to unity in the remainder of this Appendix. 
The power-law index is restricted to /3 G [— ^ , 1] with 
(3=1 and — ^ corresponding to, respectively, the har- 
monic potential and that created by a point mass at the 
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Table Bl: Relations for circular orbits in power-law poten- 
tials (B2) in units that imply _Ro = 1 = '^o- 
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origin. The fundamental quantities of circular orbits ex- 
pressed as functions of radius, angular momentum, or 
energy are listed in Table Bl, while the relations for k 
follow from 



7=^2/(1 + /?). 



(B3) 



In general, the radial action J^ cannot be expressed 
in terms of elementary functions of {E,L); however, for 
the purely radial orbits with L — 0, Jn may be given in 
closed form: 
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with b' = b/xa. Here, I have assumed that b was chosen 
according to (|2J), which gives dK/dE^ —Sk^^x^'^ and 

(CIO) 
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For the isochrone potential (|2^), g = 0, while for a flat 
rotation curve, 
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The Isochrone Approximation 



are 



The approximation for -R(t) and 0(t) resulting from the 
the isochrone approximation of § ^.2.1 
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